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Models which associate the origin of the ultra-high energy cosmic rays (UHECR) 
with decays of relic superheavy particles predict the anisotropy of UHECR flux 
toward the Galactic center. We use the existing SUGAR data, which covered the 
Galactic center best so far, to look for such a signal and limit the fraction of UHECR 
produced by this mechanism. The absence of anisotropy toward the Galactic center 
in the SUGAR data implies at 95% confidence level that the fraction of SHDM- 
related cosmic rays should be less than 20% at E > 1 x 10 19 eV and less than 50% 
at E > 4 x 10 19 eV. 



I. INTRODUCTION 

The origin of ultra-high energy cosmic rays (UHECR) is a long-standing puzzle. Observa- 
tions of cosmic rays with energy above 10 19 eV reveal that the GZK cutoff p| in the energy 
spectrum is not observed by AGASA j2] (though it seems to be consistent with HiRes jj), 
and that the arrival directions of cosmic rays are quite isotropic at large scales but cluster 
at small angles^. This is a perplexing situation for seeking UHECR sources. The isotropic 
distribution of UHECR arrival directions, in the absence of strong magnetic field which can 
produce it, either Galactic or extraGalactic, implies that the sources are distributed over a 
few hundred Mpc for which the matter distribution looks homogeneous, while the absence 
of the GZK cutoff implies that the UHECR sources are within the GZK radius which is 
about 50 Mpc. The observed correlations of UHECR with BL Lacertae objects do not 
make the situation easier: although perfectly consistent with clustering, BL Lacs as sources 
of UHECR cannot explain super-GZK events under assumption of standard physics jfj. 

A possible way to reconcile these contradictory facts is to assume that, apart from 
UHECR accelerated in BL Lacs, there exists a second, GZK-violating component orig- 
inating from decays or annihilation of superheavy particles within the GZK radius 0, 
These superheavy particles are produced in the early universe [HI and constitute a fraction 
of cold dark matter. The latter is concentrated in galactic halos. Since we are 8.5 kpc away 
from the Galactic center, the clear signature of superheavy particle decays is the signifi- 
cant excess of the UHECR flux from the direction of the Galactic center [11]. This kind 
of anisotropy toward the Galactic center (or toward any other strong source) has not been 
observed by the Northern hemisphere arrays 2, h|. However, Northern arrays do not see 



the Galactic center directly, and the small statistics of cosmic rays with energy above the 
GZK cutoff makes it difficult to draw strong constraints. 
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In this paper we use the SUGAR data |l4| to examine the viability of the superheavy dark 
matter (SHDM) scenario of UHECR, for SUGAR is the unique experiment which operated in 
the Southern hemisphere before Pierre Auger Observatory and covered the Galactic center. 
We use the Kolmogorov-Smirnov (K-S) test to check the consistency of the SUGAR data with 
the arrival direction distribution expected from superheavy particle decays in the Galactic 
halo. Our purpose is to limit a fraction of UHECR which are produced by this mechanism. 

We perform the analysis in two steps. In Sects. II and III we analyze the SUGAR data 
and derive bounds on the relative weight rj of the halo component in the total UHECR 
flux. These bounds are presented in Figures 3 and 4. They depend exclusively on the halo 
profile and statistics of SUGAR data. We then move in Sect. IV to the implications of these 
constraints for the superheavy dark matter models. The difference with the previous case 
is that SHDM contribution contains both Galactic and extragalactic parts. Their relative 
weight depends on a particular SHDM model (spectrum and maximum energy) and on 
absolute dark matter density in the halo of our Galaxy. Because of the uncertainties in 
these parameters we separate this question from the more robust analysis of sects. II and 
III. We will see, however, that in SHDM models the extragalactic contribution is generically 
subdominant in the energy range of interest, and therefore boudns of Sees. II and III apply 
to SHDM part of cosmic rays as a whole. We present our conclusions in Sect.V. 



II. GALACTIC HALO MODELS AND SUGAR DATA 

We assume that superheavy particles constitute a part of cold dark matter and their 
spatial distribution follows the dark matter density. We examine two different Galactic halo 
models. The first one is the isothermal (ISO) halo model ["La] , which is characterized by the 
density profile 

PmoW - (1) 

The second one is the Navarro- Frenk- White (NFW) model 0|, which is obtained from the 
A-body simulation and thus can be more realistic 

Here R is the distance from the Galactic center, R c is the scale of the halo core, and R s is 
the characteristic scale of the NFW model. For our galaxy, the size of R c is supposed to be 
a few kpc. The value of R s is expected to be 10-20 kpc on the basis of correlation between 
the predicted total mass and the value of R s 16[. In our analysis, we leave R c and R s as 



free parameters and use R c = 4 kpc and R s = 15 kpc as prototype values. 

For cosmic rays coming from the Galactic halo, the expansion of the universe and the 
evolution of halo can be neglected. We also neglect the effects of the Galactic magnetic 
field on the propagation of UHECR and will make a comment on it afterward. Since we 
consider the halo model which is spherically symmetric about the Galactic center and the 
distribution of sources follows the halo profile, the cosmic ray flux can be calculated as 

m) = —\ L(Jr 2 -2rRo cos 9 + J2g) dr, (3) 

47T JO 

where L(R) is the luminosity of cosmic ray sources, 9 is the angle measured from the Galactic 
center direction, which is related to the Galactic coordinates I and b by cos^ = cos b cos/, 
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r max(#) = RqCos9 + V -Rmax — -^o s i n ® ■> -^max is the extension of the Galactic halo, and 
Ro ~ 8.5 kpc is the distance to the Earth from the Galactic center. In the case of decay, 
the luminosity L(R) is proportional to the number density of superheavy particles, that is 
L(R) oc nsu(R) oc pnaio(R)- 111 the case of annihilation, L is proportional to ng H . Since 
the density profile has a peak at the Galactic center, the annihilation case predicts stronger 
anisotropy than the decay case and leads to stronger constraints on the halo contribution 
to UHECRs. In the following, we will consider the decay case only. 

To obtain the arrival direction distribution predicted from the Galactic halo model, we 
use the Monte-Carlo simulation based on the flux Q. In simulating cosmic ray data detected 
by an array of detectors, we also need to consider the exposure function of the detector array. 
Here, we use the simple geometrical exposure function which is a function of declination 5 
only 

h(5) = — [sin ct m cos A cos 5 + a m sin A sin 8) (4) 

7T 

where A is the latitude of the detector array, 9 m is the zenith angle cut, 

cos 9 m — sin A sin 5 . . 

£ — T 7 1 (5) 

cos A cos o 

and a m = for £ > 1, a m = tv for £ < — 1 and a m = cos -1 £ otherwise. 

The SUGAR array, located in the Southern hemisphere, has the latitude A = —30.53°. 
It covered the Galactic center direction best so far, because other detectors are all located 
in Northern hemisphere and cover the Galactic center marginally at best. Therefore, the 
SUGAR data is most suited at the moment for the study of cosmic rays coming from the 
Galactic halo, though in the near future Pierre Auger array will provide much better data. 
The published SUGAR data contain events with zenith angle up to 72° 111. We use the 



data with zenith angle smaller than 60° considering increasing uncertainty at large zenith 
angles. Figure 1 shows arrival directions of 316 SUGAR events with zenith angles z < 60° 
and energy E > 1.0 X 10 19 eV (upper panel), and simulated data obtained using the flux 
(|2J) for two halo models (lower panels). Note that the NFW model with R s = 15 kpc and 
the isothermal model with R c = 4 kpc predict similar concentration of events around the 
Galactic center. 



III. KOLMOGOROV-SMIRNOV TEST 

The most direct way to compare the arrival direction distribution of the real data and 
that predicted by the superheavy dark matter model is to use the two-dimensional K-S test 
over the pair (a, 5) or (/, b). Here, we use a different method which is simpler and gives better 
results. Since the cosmic rays in the halo model are distributed symmetrically around the 
Galactic center, we analyze the distribution of cosmic rays over the angle 9 between their 
arrival directions and the direction to the Galactic center. We use the one- dimensional K-S 
test to compare the SUGAR data and model predictions. Since one- dimensional K-S test is 
reparametrization-invariant, we use x = cos 9 which is more convenient in computation. 

The K-S test makes use of the cumulative probability distribution function, which is 
defined by 

S N (x) = — J2 e ( x i ~ x )i ( 6 ) 
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FIG. 1: (a) The sky map of arrival directions of 316 SUGAR events with zenith angle z < 60° and 
energy E > 1.0 x 10 19 eV. (b) and (c) show 316 x 50 simulated events for the isothermal halo model 
with R c = 4 kpc and the NFW halo model with R s = 15 kpc, respectively. Galactic coordinates 
are used; the solid line is the boundary of the SUGAR acceptance with the zenith angle cut of 60°. 



COS 



FIG. 2: Cumulative probability distribution functions used in Kolmogorov-Smirnov test of 316 
SUGAR data (z < 60° and E > 1 x 10 19 eV, solid), 4 x 316 simulated data of isotropic distribution 
(dash-dotted), the isothermal halo model (R c = 4 kpc, dashed), and the NFW halo model (R s = 15 
kpc, dotted). 

where N is the number of data points. For the distribution of arrival directions expected 
in the Galactic halo model, we generate large number M of simulated events and obtain 
corresponding cumulative distribution Sm(x). We then calculate the K-S statistic 



which determines the probability that the two sets of events are drawn from the same 
distribution |l7| . Small probabilities indicate that the two sets significantly differ. Figure 2 
shows the cumulative distribution functions used in the K-S test. The difference between 
the SUGAR data and halo models is clearly seen (it can be expected already from Figure 1). 

As mentioned in the Introduction, the SHDM model does not explain clustering and 
correlations of UHECR with BL Lacs, so another (extragalactic) component is needed which 
is isotropic at large angles. The SHDM model itself contains the isotropic extragalactic 
contribution coming from SHDM decays in other galaxies. Thus, let us consider a generic 
model where UHECR flux consists of the isotropic and Galactic halo components (their 
relative strength in the SHDM models is discussed in the next section). An important 
quantity is the fraction r) of cosmic rays coming from the Galactic halo, 



D = max|S , 7v(x) — Sm{%)\ 



(7) 



J Enii J G (E)dE 



(8) 



L (f G (E) + ME)) dE 
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FIG. 3: Probability that all cosmic rays come from the Galactic halo (77 = 1). Up: dependence on 
energy cut for different values of R c (isothermal halo) and R s (NFW halo). Down: dependence on 
halo parameter R c<s for different values of energy cut. 

where fc{E) and fu(E) are the energy spectra of cosmic rays from the Galactic halo and 
from the uniform background, respectively. We consider 77 as a function of the lower energy 
cut E min . When two spectra fc(E) and fy(E) have quite different shape, r](E min ) depends 
sensitively on the value of E min . Instead of following a specific model to evaluate r](E min ), 
we use rj as a parameter and calculate the probability to get the SUGAR data distribution. 
In this way bounds on 7/ at a given energy cut can be obtained. 

Figure 3 shows the dependence of the probability that all cosmic rays come from the 
Galactic halo (77 = 1) as a function of the energy cut .Emm and the halo parameters R c and 
R s . The dependence on R s in the case of NFW halo profile is rather weak. At realistic values 
of parameters R c ~ 4 kpc and R s ~ 15 kpc the two halo models give similar constraints. 
Dominant fraction of halo component is excluded up to energies as high as 5 x 10 19 eV. 
For higher energies, the constraints become progressively weaker; the reason may be much 
smaller statistics. 

Figure 4 shows the lines of constant probability (corresponding to la - 4a levels) in the 
plane of the model parameters R c , R s and the fraction of halo component 77. Two energy 
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FIG. 4: Exclusion plots (lines of constant probability) for the two-component model described in 
the text on the r\ - R c (isothermal halo, left panels) and rj - R s (NFW halo, right panels) plane. 
Two energy cuts are presented, E m [ n = 1 x 10 19 eV (upper panels) and E m i n = 4 x 10 19 eV (lower 
panels). Symbols •, +, o and • represent la - 4<r lines, respectively. 



cuts are represented: E min = 1 x 10 19 eV (upper panels) and E min = 4 x 10 19 eV (lower 
panels). In the case of NFW halo profile, isotropy of SUGAR data at E > 1 x 10 19 eV 
excludes fractions of halo component larger than 15% at 2a level, while for E > 4 x 10 19 eV 
the halo fraction has to be less than 45%. Constraints are similar for isothermal halo profile 
eq.(JD) at R c ~ 4 kpc. 

Finally, let us comment on how robust are the results presented here. One thing one 
should check is the influence of the Galactic magnetic field during the propagation. This can 
be important because significant portion of cosmic rays from the Galactic halo are produced 
near the Galactic center and experience strong magnetic fields. To see the effect of this 
spread on our results, we added a Gaussian random deflections with mean deviation up to 
10° in the simulated data set. This effect also mimics possible poor angular resolution of 
SUGAR events. We found that results presented above are stable against this modification. 
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The maximum change occurs at E min « 1 x 10 19 eV and r\ « 1. For the moderate energy 
cut or with the uniform background component added, its effect becomes less than 10% and 
can be completely neglected. 



IV. APPLICATION TO SHDM MODEL 

In order to see what the above constraints imply for a given SHDM model, one has 
to estimate the strength of the halo (anisotropic) component relative to the extragalactic 
(isotropic) one. The latter consists of a part related to SHDM decays outside of the Galaxy 
and a part due to astrophysical sources. Let us estimate the former. The Galactic and 
extragalactic contributions in the SHDM model can be written as 



F G {E) = C J d 3 x 



n(x) 



(R-x) 



F cxt (E) = C J d 3 x^w(r). 



Here n and n(x) are the average and halo dark matter densities, R is the vector from 
Galactic center to the Earth, while w(r) is the weight function which accounts for the flux 
attenuation at large distances. The proportionality coefficient C which accounts for the 
efficiency of UHECR production from SHDM is the same in both cases. It therefore cancels 
in the ratio of the two contributions. 

The extragalactic part can be written in the form 



ext 



2nCnR e fi, 



where by definition R e d = J drw(r). This quantity represents the size of the region from 
which the extragalactic flux is collected (for instance, for very high-energy protons R e s is 
the GZK distance). It can be estimated as the attenuation length divided by the spectral 
index. The factor 2n in this equation takes into account the fact that SUGAR sees roughly 
half of the sky. 

To estimate the halo contribution we take the NFW profile at R s = 15 kpc and calculate 
the flux coming from half of a sphere overlooking the Galactic center. The result can be 
written as follows, 

Fq ~ 2nCn x 5000 Mpc. 

When obtaining this number the ratio of the local to average dark matter density was taken 
to be 2 x 10 5 . This number is uncertain by a factor 2 or so. 
The ratio of the two components is 



ext 



ReS 



F G 5000 Mpc 

The UHECR flux in the SHDM models is dominated by photons which have attenuation 
length less than 100 Mpc in the range of energies 10 19 — 10 20 eV even at very weak assump- 
tions about radio photon background |l8|. We therefore arrive at the conclusion that in the 
SHDM models the halo component dominates already at energies around 10 19 eV, i.e., quite 
below the GZK cutoff. Since the halo component dominates, the bounds on the fraction r) 
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in Figures 3 and 4 essentially apply to the total fraction of UHECR which originate from 
decays of SHDM particles. 

The situation becomes even better at high energies where the extragalactic component 
dies away and the halo contribution reaches 100%. In this energy range ~ 60 uniformly 
distributed events would be enough to rule out the SHDM model at 4 — 5a level following 
the technique developed here. 

V. CONCLUSIONS 

We have used the SUGAR data to check the viability of the hypothesis that a part of the 
highest-energy cosmic rays originate from decays of relic SHDM particles. We found that 
the absence of the anisotropy toward the Galactic center in the SUGAR data imposes strong 
constraints on corresponding fraction of UHECR. At 95% confidence level, the fraction of 
SHDM-related cosmic rays should be less than 20% and 50% at E > 1 x 10 19 eV and 
E > 4 x 10 19 eV, respectively. This implies that the SHDM component of UHECR, if exists, 
must have harder spectrum than the uniform one and should take over at energies higher 
than 4 x 10 19 eV. 
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